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Abstract 

We investigate baryon transport in relativistic heavy-ion collisions at energies reached at the 
CERN Super Proton Synchrotron (SPS), BNL Relativistic Heavy-Ion Collider (RHIC), and 
CERN's Large Hadron Collider (LHC) in the model of saturation. An analytical scaling law 
is derived within the color glass condensate framework based on small-coupling QCD. Transverse 
momentum spectra, net-baryon rapidity distributions, and their energy, mass, and centrality de- 
pendencies are well described. In comparison with RHIC data in Au + Au collisions at ^snn = 
62.4 and 200 GeV, the gradual approach to the gluon saturation regime is investigated and limits 
for the saturation-scale exponent are determined. Predictions for net-baryon rapidity spectra and 
the mean rapidity loss in central Pb + Pb collisions at LHC energies of ^/snn = 5.52 TeV are 
made. 

PACS numbers: 24.85. +p, 25.75.Dw, 12.38.Mh 
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I. INTRODUCTION 



Gluon saturation has been the focal point of important and interesting particle-physics 
investigations for many years. Its observation would allow to access a new regime of quantum 
chromodynamics where high-density gluons form a coherent state. In regions of large parton 
densities the physics is governed by a single hard scale, Q s ^> Aqcd, which increases with 
energy and thus allows the use of small-coupling techniques l). In this regime, gluon 
recombination starts to compete with the exponentially increasing gluon splitting and the 
gluon distribution function is expected to saturate. 

At the HERA, some evidence for gluon saturation in the proton was found in deep inelastic 
e + p collisions at high energy and low values of Bjorken-x, but the results are still open to 
interpretation The existence of geometric scaling as predicted by the color glass theory 
as an approach to saturation physics was indeed confirmed, constituting the most important 
evidence for saturation so far 3j. 

Since the saturation scale is enhanced by a factor A 1 / 3 in heavy ions, as compared to 
protons, it is natural to investigate saturation in relativistic heavy-ion collisions, as has 
been done by many authors 4|. Here theoretical QCD-based approaches have, so far, usually 
concentrated on charged-hadron production and, in the central rapidity region, a reasonable 

nnfiQ 

understanding was achieved in the color glass condensate framework [1|, Ha, la, 17[] through 
inclusive gluon production p, [9|. 

In p + A collisions the nuclear wave function is being probed at small x in the forward 
direction and it should show saturation below a characteristic value of x. Hence, the effect of 
gluon condensation and quantum evolution should be measurable, as was attempted by the 
3RAHMS Collaboration at the Relativistic Heavy- Ion Collider (RHIC) in d + Au collisions 
lo| . Measurements of the nuclear modification factor showed the Cronin enhancement 
at midrapidity, but the modification factor is suppressed at forward rapidities. This is in 



12|. 



qualitative agreement with the color-glass predictions 11, 

Experimental heavy-ion investigations at the Large Hadron Collider (LHC) concentrate 
on the midrapidity region since ALICE [3] covers rapidities up to \y\ = 2. It provides 
measurements of lower x values than before — down to 10~ 5 — at an energy scale that is high 
enough to provide crucial tests of gluon saturation. 

In the present phenomenological investigation we use the transverse momentum spectra 
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and rapidity distributions of net baryons (B — B) in relativistic heavy-ion collisions 
as a testing ground for saturation physics. A primary account of our approach was given in 
Ref. [if]]. It made use of the color-glass framework []], 0, [fl Q, [2, [9}. Related approaches of 



other authors to the net-baryon problem with respect to gluon saturation are Refs. 



13, 
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We compare with data that were obtained by scaling distributions of identified net protons at 



Super Proton Synchrotron (SPS) and RHIC energies, and extrapolate to LHC energies. 
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problem was treated with different approaches based on QCD without saturation 

In A + A collisions, two distinct and symmetric peaks with respect to rapidity y occur 
at SPS energies [22| and beyond. The rapidity separation between the peaks increases with 
energy and decreases with increasing mass number, A, reflecting larger baryon stopping 
for heavier nuclei, as was investigated phenomenologically in the nonequilibrium-statistical 



relativistic diffusion model 



23 



241 ] . In this work we show how the evolution of the peaks 



can be linked to saturation physics. 

The net-baryon number is essentially transported by valence quarks that probe the sat- 
uration regime in the target by multiple scatterings. During the collision, the fast valence 
quarks in one nucleus scatter in the other nucleus by exchanging soft gluons, leading to 
their redistribution in rapidity space. We take advantage of the fact that the valence quark 
parton distribution is well known at large x, which corresponds to the forward and backward 
rapidity region, to access the gluon distribution at small x in the target nucleus. Therefore, 
this picture provides a clean probe of the unintegrated gluon distribution, {p(x,pr), at small 
x in the saturation regime. Here pt is the transverse momentum transfer. 

In particular, we use net-baryon rapidity distributions in central relativistic heavy-ion 
collisions from SPS to LHC energies to probe saturation physics through their energy and 
mass-number dependence on a geometric scaling variable r. 

We discuss the net-baryon rapidity distributions from SPS to LHC energies in Sec. II. 
Their dependence on the scaling variable r is investigated and the position of the fragmen- 
tation peak as a function of the saturation-scale exponent A is derived. Possible conclusions 
about saturation from the midrapidity valley in net baryons are drawn and the effect of 
quark fragmentation into hadrons is considered. 

In Sec. Ill, theoretical results are compared with data on transverse momentum spectra 
for Au + Au at RHIC energies of 62.4 and 200 GeV and with rapidity spectra at SPS energies 
of 17.3 GeV, as well as at RHIC energies. Predictions for the rapidity distribution of net 
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baryons at LHC energies are made. Net-kaon rapidity distributions are also discussed. 
Centrality and system-size dependence are investigated, and the mean rapidity loss as a 
function of beam rapidity, or energy, is discussed. Finally, conclusions are drawn in Sec. IV. 



II. RAPIDITY DISTRIBUTIONS IN THE COLOR GLASS CONDENSATE 

The differential cross section for valence quark production in a high-energy nucleus- 



nucleus collision reads 
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d^y = (2^ Xiq ^ Q ^ V (X2 ' PT) ' (1) 
where pt is the transverse momentum of the produced quark and y its rapidity. The lon- 
gitudinal momentum fractions carried, respectively, by the valence quark in the projectile 
and the soft gluon in the target are x\ = pT/\^sexp(y) and x<i = pr/ y/sexp(— y). The 
contribution of valence quarks in the other beam nucleus is added incoherently by changing 
y — > —y. The valence quark distribution of a nucleus, q v = q — q, is assumed to be given by 
the sum of valence quark distributions q Vt N of individual nucleons, q v = Aq Vt ^, where A is 
the atomic mass number. 

The factorization scale is usually set equal to the transverse momentum, Qf = Pt- Since 
the valence quark parton distribution is weakly dependent on Qf we neglect it in the following 
discussion. The gluon distribution is related to the forward dipole scattering amplitude 



J\f(x,rT) 291 ] . for a quark dipole of transverse size rx, through the Fourier transform 

(f(x,p T ) = 27ip T J r T dr T J\f(x,r T )J (r T p T ). (2) 

Assuming the rapidity distribution for net-baryons is proportional to the valence quark 
rapidity distribution up to a constant factor of C, we obtain, by integrating over pt, 

dN C f d 2 p T 

^ = J^J -jrxM*i,Qf)<P(x2,Pr). (3) 

We show that this is indeed a good approximation, at high energy, in the section devoted 
to the issue of fragmentation functions. 

A. Scaling property of the net-baryon distribution 

Since the valence quark distribution is well known, net-baryon production represents a 
good observable to probe saturation effects in high-energy heavy-ion collisions through the 
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FIG. 1: Projectile contributions to the net-baryon rapidity distribution plotted as a function of the 
variable y' = y — yt (upper panel) for central Au + Au collisions at RHIC energies of ■ v /s/vtv = 62.4 
GeV (dotted curve) and 200 GeV (dashed curve) and for Pb + Pb at LHC energies of ^/snn = 
5520 GeV (solid curve). The limiting fragmentation property is obviously broken for net baryons. 
The lower panel shows the equivalent curves for our variable y' = y — y\ ) /{l + A) with A = 0.2; here 
scaling is fulfilled in the peak region. 

respective gluon distributions. One important prediction of the color glass condensate theory 
is geometric scaling: The gluon distribution depends on x and Pt only through the scaling 
variable p^/Q 2 s {x), where Q 2 s (x) = A l /^Ql x~ x ; the scale Qo sets the dimension. Geometric 



scaling was confirmed experimentally at HERA [3]. The fit value A = 0.2-0.3 agrees with 



theoretica 
results 30 



estimates based on next-to-leading order Balitskii-Fadin-Kuraev-Lipatov (BFKL) 



31] . To show that, in the high-energy limit, the net-baryon distribution reflects 



the geometric scaling of the gluon distribution, we perform the following change of variables: 



— 2v 2 2 — 2v 

x = Xi, x 2 = x e y , p T = x s e y . 



(4) 
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Thus, we rewrite Eq. ([3} as 



dN , , C f 1 dx 



T 



xq v (x) <^(x 2 +V), (5) 

where r = \n(s/Ql) — In A 1 / 3 — 2(1 + A) y is the corresponding scaling variable. Hence, the 
net-baryon multiplicity is only a function of a single scaling variable r, which relates the 
energy dependence to the rapidity and mass-number dependence. From the equation for the 
isolines, r = const, we obtain the evolution of the position of the fragmentation peak. This 
reflects the interplay of the valence quark distribution, peaked at x\ ~ 0.2, with the gluon 
distribution peaked at pt ~ Q s , in the forward region with respect to the variables of the 
problem 

2/peak = "j ~ r {yb ~ In A 1/6 ) + const, (6) 
1 + A 

where yb = (1/2) \n[(E + pi)/(E — pi)] — ln(i/i/m) is the beam rapidity at beam energy E 
and longitudinal momentum p^ with the nucleon mass m. 

In Fig. [U we show the peak positions as given by Eq. ([6]) for different incident energies. Fl 
Hence, different values of yb, as functions of y' with y' = y — yb in the upper panel, and 
y' = y — y&/(l + A) in the lower panel, where y' is linearly related to our scaling variable 
r through r = —2(1 + \)y' are shown. In this calculation we include baryon mass effects 
by the replacement pt — > \/pt + m 2 . We then observe clear violations of the y — yb scaling 
in the upper panel. Hence, according to our model, limiting fragmentation phenomena as 



observed in particle production [33|], where the rapidity distribution scales as y — yb, is 
violated in net-baryon rapidity distributions, which should instead exhibit a scaling with 
y' = y — yb/{l + A)This is shown in the lower panel, where scaling is fulfilled in the peak 
region. The deviations outside the peak region are largely due to mass effects. To draw this 
conclusion we do not yet need to specify the form of the gluon distribution, although for 
numerical calculations we must specify one. The details of the computations are given in 

Sec. dm 



The compilation of the SPS data at ^/s =17.3 GeV |22j and RHIC at = 62.4 and 200 
GeV [l5| provides an opportunity to verify the scaling law in Eq. (J6j). The scaling properties 
of net-baryon rapidity distributions were investigated by the BRAHMS Collaboration for 
the one-nucleus contribution (see Ref. [15] for more details about the subtraction procedure 
of the projectile contribution). Despite the lack of data in the peak region at RHIC, the data 
seem to exhibit a scaling with y — yb compatible with the limiting fragmentation picture. 
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A new and more refined analysis of SPS data was recently done by the NA49 Collaboration 



271 ]. showing a slight discrepancy with the previous data 22J. The new data are shown as 



black squares in Fig. [2]for y' = y—yb (upper panel) and for the scaling with y' = y—yb/ (1+A) F2 
(lower panel) with A = 0.2. The present data do not allow us to distinguish between the two 
different scaling laws. The value of the saturation-scale exponent A = 0.2 was determined 
by recent calculations [18j in the saturation picture, including running coupling effects, 
in agreement with particle production at RHIC. To improve the agreement with the new 
SPS data as compared to the analysis reported in Ref. |l6[], which was based on the old 
NA49 data, we increased A from 0.15 to 0.2 in the present work. Note, however, that 
the applicability of the color glass condensate (CGC) picture may be questionable at the 
relatively low SPS energies. 

In the peak region, the average x in the projectile is x ~ 0.2-0.3, which corresponds to the 
average momentum fraction carried by a valence quark. In the target, x = (0.2-0.3) e~ 2y p cak , 
it decreases with increasing energy. In this kinematic regime we have a natural intrinsic 
hard momentum, the saturation scale Q s . This justifies the use of small-coupling techniques 
in QCD for calculating integrated yields 38( . The effects of the medium are expected to be 
small at forward rapidity since the fast moving valence quarks escape the interaction zone 
quickly. A detailed measurement of the peak region will then enable us to reconstruct the 
gluon distribution from Eq. ([3]). 



B. Midrapidity valley in net protons 

It is worthwhile to analytically investigate some limits of Eq. ([3]). In this section, we 
derive a parametric formula for the region of small Xi (away from the peak), corresponding 
to the midrapidity valley (y ~ 0). In this kinematic regime the valence quark distribution 
behaves where A ~ 0.5 is the intercept of the Regge trajectory [17] . which 

allows us to perform analytic calculations. 

First, let us recall that the unintegrated gluon distribution is peaked at qx = Q s , or 
X\ = exp (— t/2 + A), reflecting the fact that most of the gluons sit at this value. Therefore, 
we expect dN/dy ~ xiq((xi)), with (xi) = (Q s ) / \fs exp(y). Recalling that Q 2 S = A 1//3 QqX^ a 
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FIG. 2: Measured projectile contributions to the net-baryon rapidity distribution plotted as a 
function of y' = y — y& (upper panel) and y' = y — y^j '(1 + A) (lower panel), with a saturation-scale 
exponent A = 0.2. Black squares are preliminary new Pb + Pb SPS data at ^snn = 17.3 GeV 



, open squares are old Pb + Pb NA49 data 



22], circles are Au + Au BRAHMS data at 62.4 



ii. 



GeV [HI, and triangles are BRAHMS data at 200 GeV 
and X2 = x\ exp(— 2y), we can solve the equation for (xi), yielding 

(*i> 







[9 1 + A 1 




) exp 


[ 2 2 + X y \ 



Finally, we obtain 



1 dN (A X ' 6 Q { x A/{1+|) 



oc 



cosh 



2 + A y 



(7) 



(8) 



A dy " \ ^/s 

This property reflects the following asymptotic behavior at large r, and x <C 1, of the 
gluon distribution: 

lim <p(x, t) oc S(x 2+x e T - 1). (9) 

r— »oo 

This approximation is valid as long as r is large (and iCl), typically in the midrapidity 
valley at asymptotic energies. For A = we recover Eq. (80) in Ref. 17] plus the A 
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dependence: 

1 dN fA l l & \ A . . 

oc — =■ cosh(2A?/). (10) 



A dy \ 
C. The fragmentation function 

In this section we investigate the effect of quark fragmentation on our conclusions about 
geometric scaling. By including the fragmentation function of a quark into net baryons, 
D(z) = D/± B / q (z) = D B / q (z) — D B / q (z) [z = E B /E q being the fraction of the quark en- 
ergy carried by the baryon fragment) and using isospin symmetry D(z) = D&B/q{z) = 
—D& B /q(z), the cross section for the production of a hadron of transverse momentum pt at 
rapidity y reads 

where = \pPr + m2 l ' z is the quark momentum, xp = \/p\ + m 2 / y/s exp(y) is the 
Feynman-x, and Xi = qx / 'y/s exp(y) , x 2 = qr/ \/sexp(—y). Integrating Eq. ffTTT) over p T up 
to the kinematic boundary p max = y/s e~ y imposed by xf < 1, and inverting the order of 
the integrals, yields 
dN 



dy Jo (2tt; 



2 



>-F 



(27T) 













Qt 


/"9max 


d 2 q T 


) / 

J m/z 


<1t 



/l Mmax d?qT 

dzD(z) / xi<?„(xi) if {x 2 , qr) , 

.0 Jm/z Qt 

(12) 



with z = m/y / iexp(y) and g max = y/sexp(—y). Now, if we assume geometric scaling we 
have 

**L = J_ f dzD(z) f - xq v (x) <p (* 2 +V) . (13) 
dy 2tt J zo J Zq/z x 

Obviously, since z depends on the energy and on the rapidity, it violates explicitly geometric 
scaling. However, in the high-energy limit, when s — > oo, or more precisely when m <C (pt), 
the lower bound of the integral can be set to and one recovers the geometric scaling formula 
of Eq. © with C = /J dzD(z). 

At this stage, we anticipate the discussion on the applicability of fragmentation functions 
for the observables of interest. The main contribution to the rapidity distribution comes 
from baryons of transverse momentum p? ~ 1 GeV, which is low enough to render the use of 
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the fragmentation functions questionable. Moreover, it was pointed out clearly by Bass and 
co-workers j^] that, for pt < 5 GeV, parton recombination dominates the hadronization 
process in particle production. Hence, the hadrons will be produced by partons of smaller 
energy instead of fragmentation of partons of larger energy as required by the fragmentation 
picture. Therefore, to simplify the discussion, we assume D(z) oc 5(1 — z) to account roughly 
for the competition between recombination and fragmentation. 

To gain more insight into the hadronization process of the valence quarks, let us recall 
the relationship between the peak position and the parton transverse momentum. At the 
peak position, the peak of the valence quark distribution x% ~ 0.2 is reflected. Hence, one 
can extract the mean transverse momentum as 

(q T ) ~ 0.2v/sexp(-?/ pcak ). (14) 

Moreover, the mean hadron transverse momentum (p?) can be extracted from the hadron 
spectra. Therefore, one obtains a measure of the energy that flows from the valence quark 
to the net protons (baryons) as 

«-S- (15) 

For (z) < 1, the hadronization process is dominated by fragmentation, in other words, 
the valence quarks lose energy by radiating gluons and, therefore, the produced baryon 
carries a fraction of the quark momentum. Furthermore, for (z) > 1, the mean momentum 
of the hadron is larger than that of the valence parton, which can be achieved by parton 



recombination 25]. 



III. THEORY VERSUS DATA 



To take into account saturation effects in the target we choose the Golec-Biernat-Wiisthoff 



model |34j for the forward dipole scattering amplitude J\f, leading to (cf. Eq. (j2j) and Ref. 



Pt ( Pt 



p{x,p T ) = In— — exp —77^ , (16) 

Qs(x) V Qs{ x )J 

in the fundamental representation of SU(3). This Gaussian form actually reflects the multi- 
ple scatterings performed by the valence quarks in the color glass. It is interesting to make 
the connection to the fitting procedure performed by the BRAHMS Collaboration when in- 
tegrating the spectra. Indeed, while the usual fitting functions used for particle production 
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are exponential in pr or have a Boltzmann form, the BRAHMS Collaboration noticed that, 
for net-proton spectra, the best fits are obtained with Gaussians in px- This corroborates 
our picture. 

First, we investigate transverse momentum distributions for net protons in comparison 
with BRAHMS data taken at different rapidities to fix Q . The spectra constrain our model 
that contains only two parameters, Qo, and the overall normalization C. This investigation 
also provides some hints regarding the transition from a coherent gluonic state (the color 
glass condensate) at low transverse momenta to incoherent partonic (jetlike) interactions at 
high transverse momenta. 

The valence quark parton distribution function (pdf ) of the nucleus is taken to be equal to 
the valence quark pdf in a nucleon times the number of participants in the nucleus. Here, we 
focus on the forward rapidity region and interpolate to midrapidity where small- x quarks 



are dominant, by matching the leading-order distributions 35_J and the Regge trajectory, 
xq v oc x - 5 , at x = 0.01 [171 ]. 

We compare the data to two calculations corresponding to different assumptions concern- 
ing hadronization of the valence quarks. First, we use fragmentation functions for valence 
quarks to net-protons 40J (the only ones that are presently available), 

D p ^{z)=N z a {l-z)\ (17) 

with N = 520142, a = 11.6, and b = 6.74 fitted at Q = 1.4 GeV. We can neglect QCD 
evolution in the range of interest. This is represented by dashed curves in Figs. |3] and HI F3,4 
We refer to this result as the F-model (for fragmentation model). In the NF-model (no- 
fragmentation), we assume that, on average, the energy of the parton is equal to that of 
the produced baryon, namely D(z) oc 5(1 — z). This assumption is based on the fact that 



hadronization at low px is poorly known and it was proven 25] that it is dominated by 
recombination of partons instead of fragmentation, which will imply z > 1. 

Hence, there should be a competition between the two phenomena, fragmentation and re- 
combination. Namely, due to the high parton density characteristic of high-energy heavy-ion 
collisions, the parton showers described in the fragmentation picture can overlap. Hence, the 
recombination of partons from different showers can occur in addition to the recombination 



of shower partons with thermal partons 25|, |26 |. 



To fix the scale Q we use the forward-rapidity net-proton spectra at 200 GeV where we 
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expect our model to be valid: high energy and large rapidity. Whereas Qo is fixed once for 
all the available energies, the normalization is tuned at each energy to fit the data owing to 
our lack of knowledge on the hadronization process. Moreover, our calculation is at leading 
order, therefore, a K-factor is required. 

For A = 0.2, we fix Q 2 = 0.04 GeV 2 in the NF-model (0.1 GeV 2 in the F-model), leading 
to Ql = 0.6 GeV 2 (1.5 GeV 2 in the F-model) at x = 0.01. When comparing to investigations 
of charged-hadron production, such as in Refs. j^l, 12], which involve the gluon distribution 
in the adjoint representation of SU(3), we must consider a rescaling of our net-baryon Q 2 S 
by the color factor Nq/Cf with Cf = (Nq — l)/2Nc and Nc = 3, corresponding to a factor 
9/4. 

A. Transverse momentum spectra and rapidity distributions at 200 and 62.4 GeV 

In Figs. [3] and HJ the pt spectra at 200 and 62.4 GeV are shown for a large rapidity 
range. Results for the NF-model (full curves) are compared with the F-model (dashed 
curves). The increasing discrepancy between the data and both models toward midrapidity 
is expected because of the reduction of the window for saturation effects with decreasing 
energy and rapidity. In addition, the quark-nucleus cross section is valid in the eikonal 
approximation: It is meant to be applicable at forward rapidity where the valence-quark 
longitudinal momentum fraction is X\ ~ 1. Another source of uncertainty can result from 
an additional rescattering of the valence quark in the quark-gluon plasma leading to an 
enhancement of the mean broadening, essentially at midrapidity, where medium effects 
are expected to be important, 

(Pt) ^Qs^ (Pt) ^Qs + (PT)med, (18) 

resulting in a harder spectrum at midrapidity and better agreement with the data. 

The integrated net-proton rapidity distributions, scaled by a factor of 2.05 [jj]] to ob- 
tain the net-baryon distributions, are shown in Fig. [51 Whereas the NF-model (full line) F5 
describes the data well, the implementation of the fragmentation function gives an unsatis- 
factory result. The estimated numbers of participants are 390, 315, and 357 for -y/s = 17.3, 



62.4, and 200 GeV, respectively 15]. From the integration of the rapidity distribution of the 



NF-model about 15% of the participant baryons are missing, most of them in the tails of the 
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FIG. 3: Net-proton spectra for Au + Au at RHIC energies of ^/snn = 200 GeV for different 
rapidities. From left to right, y = 0,0.9, 1.9, and 2.9. The triangles represent BRAHMS data 



20 | . The saturation-scale 



141]: at y = 2.9, more recent preliminary BRAHMS data are shown 
exponent in the calculations is A = 0.2. Full lines correspond to the no-fragmentation hypothesis 
with Qq = 0.04 GeV 2 , whereas the dashed lines include fragmentation with Qq = 0.1 GeV 2 . 

distributions, which are too steep in our high-energy model. Hence, the NF-model accounts 
for about 85% of the estimated baryon number, whereas the F-model, which includes frag- 
mentation, accounts for significantly fewer baryons. We believe that a better description of 
the hadronization process, in particular by incorporating recombination, will improve the 
overall agreement with the data. 

We extrapolate our results for the NF-model to central (0%-5%) Pb + Pb collisions 
at LHC energies of 5.52 TeV in Fig. [6j where the theoretical distribution is shown for F6 
A = 0.2 and = 0.04 GeV 2 . Since ALICE measures identified protons and antiprotons 
only in the rapidity range \y\ < 2, a direct observation of the fragmentation peaks, and of 

;he dependence of the peak position on the saturation-scale exponent as discussed in Ref. 

1^ | . will depend on future extensions of the LHC heavy-ion detectors to enhance forward 
capabilities. 

A further test of our model, and of the value of the deduced saturation scale, is provided 
by comparing it with net-kaon distribution functions and, in particular, dN/ dy for net kaons. 
The bulk of the produced K + and K~ mesons are due to inclusive gluon production in the 
midrapidity source. The net kaon distribution ( K + — K~), however, is essentially due to 



13 (November 23, 2009) 



10 



> 

£ io 0[ 

73 -2 " 
73 



1 1 1 1 





y=0.65 






vs 








i 











1 2 








1 2 1 2 1 2 3 



P T [GeV] 



FIG. 4: Net-proton spectra for Au + Au at RHIC energies of y / s/v r /v = 62.4 GeV for different 
rapidities. From left to right, y = 0, 0.65, 2.3, and 3. The squares represent the BRAHMS data 
151 ] . The full lines correspond to the no-fragmentation hypothesis with Qq = 0.04 GeV 2 and the 
dashed lines include fragmentation with Qq = 0.1 GeV 2 . The saturation-scale exponent is A = 0.2 
in all cases. 



the interaction of fast valence quarks with soft gluons in the target, just as in the case of 
the net-proton distribution that we have discussed thus far. 

The net-kaon distribution that is obtained from the 200-GeV BRAHMS data is shown in 
Fig. [7J The data show the same trend as the net-baryon results. For kaons, the net-charge F7 
content is 44. (In the total charge balance of the collision, we must consider all particle 
species. Since negative pions are slightly more abundant than positive pions, the integral of 
net-proton and kaon distributions may actually exceed 158.) The calculation, as shown in 
Fig. [7J agrees well with these data. 



B. Centrality and system-size dependence 

We also investigate the centrality dependence of the net-baryon distribution. Formally, in 
Sec. Ull we show the dependence of the rapidity distribution on the mass number A through 
the saturation scale Q s oc A 1 ^. The centrality dependence of particle production is essen- 
tially determined by the number of participants and hence we now make the replacement 
^4 — > N part . 
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FIG. 5: Rapidity distribution of net baryons in central (0%-5%) Pb + Pb collisions at SPS energies 
°f y/sNN = 17.3 GeV (left frame). The theoretical calculation for A = 0.2 and Qq = 0.04 GeV 2 is 



compared with NA49 results that were extrapolated from the net 



proton data (open circles 



m). 



27|]. At RHIC energies of ^/sTTTv 



Black diamonds are more recent preliminary NA49 data points 
62.4 GeV (middle frame, 0%-10%) and 200 GeV (right frame, 0%-5%) for central Au + Au, our 
corresponding theoretical results are compared with BRAHMS net baryon data (circles) [lj[ 151. 
At 200 GeV, triangles are preliminary scaled BRAHMS net proton data points for 0%— 10% 32J. 
The full lines correspond to the no- fragmentation hypothesis with Qq = 0.04 GeV 2 and the dashed 
lines include fragmentation with Qq = 0.1 GeV 2 . Arrows indicate the beam rapidities. 



100^ Pb+Pb 
LHC 



5.52 TeV 




FIG. 6: Rapidity distribution of net baryons in 0%— 5% central Pb + Pb collisions at LHC energies 
of y7F/v7v = 5.52 TeV. The theoretical distribution is shown for A = 0.2 and Qq = 0.04 GeV 2 . 



Let us recall first a general feature of net-baryon distributions as a function of system 
size A: More stopping is observed for larger systems. In other words, it corresponds to the 
shift of the fragmentation peaks toward midrapidity with increasing A. To quantify this, we 
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FIG. 7: Net-kaon rapidity distribution as obtained from our model in comparison with the 
BRAHMS data Q] for central (0%-5%) Au + Au at ^/saw = 200 GeV. The calculation is for 
A = 0.2 and Q§ = 0.04 GeV 2 . 



consider two experiments involving different nuclei, A\ and A 2 such that A\ > A% 3> 1, then 
according to Eq. ([6]) for the net-baryon peak position, we obtain for the rapidity difference 
of the peaks for two different mass numbers: 



A|/peak = 2/peak(A) ~ 2/peak(^2 



1 + A \A 2i 



(19) 



where the negative sign reflects the increasing stopping with increasing A. Let us consider a 
central gold-gold collision with iV part — 350 = A\ and a peripheral one with iVp art ~ 50 = A 2 . 
Using Eq. (|T9|) . for A = 0.2, we obtain for the rapidity shift Ay pcak ~ 0.27. Hence, provided 
the measurements of the peak region are precise, we access the value of A by measuring the 
peak shift. Unfortunately, the data are not completely conclusive yet. In Fig. [SJ we show F8,9 
that the data for Pb + Pb collisions at 17.3 GeV with N paxt = 352 in the measured rapidity 
range (390 if one extrapolates to full rapidity range) indeed scales with the data for S + S 
collisions with N paTt = 52 in the measured region [4l| according to Eq. 

We show in Fig. [9] the computation resulting from the NF-model, the centrality depen- 
dence of the rapidity distribution at 17.3 and 200 GeV. Whereas the 17.3-GeV and the more 
peripheral 200-GeV theoretical curves agree with the data, the more central 200-GeV curve 
does not reproduce the absolute magnitude of the data at forward rapidity (see also Fig. [5]). 
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FIG. 8: Pb + Pb net-baryon rapidity distribution at ^/snn = 17.3 GeV (open squares are NA49 
data [22]) and S + S rapidity distribution at ^/snn = 19.4 GeV (black diamonds are NA35 data 

1 /6 

411 ] ) plotted as functions of the scaling variable y' = y — (yt — m -^part)/(l + ^) with N paxt = 352 
for Pb + Pb and iV part = 52 for S + S. 



C. Mean rapidity loss 



With increasing energy; 



waves in rapidity space 



the peaks move apart, and the solutions behave like traveling 



361 ]. which can be probed experimentally at distinct values of the 
beam energy or the corresponding beam rapidity. We derive the peak position as a function 
of the beam rapidity as y pea ,k = vyb + const with the peak velocity v = 1/(1 + A) [cf. 
Eq. (|6])]. The position of the peak in rapidity space as a function of the beam energy 
can, in principle, be determined experimentally or at least estimated (RHIC). Theoretically, 
its evolution with energy provides a measure of the saturation-scale exponent A. Hence, a 
precise determination of the net-baryon fragmentation peak position as a function of beam 
energy will provide detailed information about the gluon saturation scale. 

The mean rapidity loss (Sy) = yj, — (y) is shown in Fig. [TDJ Our result is in agreement F10 
with the experimental values of baryon stopping that was obtained at Alternating Gradient 



Synchrotron (AGS) and SPS energies 22|, |37J. Assuming that the mean rapidity evolves 



similarly to the peak position, (y) = y pea k+ const, we show that the mean rapidity loss 
increases linearly at large y^. 



(Sy) 



A 



1 + A 



?/6 + const, 



(20) 
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FIG. 9: Rapidity distribution of net baryons in Pb + Pb collisions at SPS energies of y / s/v r /v = 17.3 



GeV (upper frame). The theoretical calculations for A = 0.2 and Qq = 0.04 GeV 2 are compared 



with recent preliminary NA49 results [27( for centralities of 0%-5 %, 12.5%-23.5%, and 33.5%- 
43.5% (top to bottom; net-proton data scaled to net baryons). At RHIC energies of ^snn = 200 
GeV (bottom frame) for Au + Au, our corresponding theoretical results without fragmentation 
are compared with preliminary BRAHMS net baryon data for centralities of 0%-10%, 10%-20%, 
20%-40%, and 40%-60% [32]. 



where the slope is related to A. Hence, the mean rapidity loss that accompanies the energy 
loss in the course of the slowing down of baryons provides a measure for A and thus a test 
for saturation physics. The gray band in Fig. [10] reflects the uncertainty of how to place the 
remaining baryons that are missing in our model. We refer to them by AN, with N = Ap art . 
The upper limit corresponds to the case where the missing baryons sit at the mean rapidity, 
roughly about the peak rapidity. Then the corrected mean rapidity loss is equivalent to the 
theoretical one, (Sy) corr = (6y). The lower limit corresponds to the case where they sit at 
the beam rapidity, (Sy) corr = (1 — AN/N)(5y). The full line is the mean value of the two 
calculations and it is in reasonable agreement with the upper limit of the data given by 
BRAHMS. 



18 (November 23, 2009) 



6 




FIG. 10: (Color online). The mean rapidity loss (Sy) as obtained from our theoretical results is 
plotted as a function of beam rapidity y\, (solid curve). The star at j/beam = 8.68 is our prediction 
for central Pb + Pb at LHC energies of y^jviV = 5.52 TeV with A = 0.2 and Q\ = 0.04 GeV 2 . 
Analysis results from AGS Au + Au data (E917, E802/E866, triangles) |37], SPS Pb + Pb data 
(NA49, square) and RHIC Au + Au data (BRAHMS, dots, with triangles as lower and upper 
limits) [p| are compared with the calculations. 

IV. CONCLUSION 

In summary, we present a saturation model for net-baryon distributions that successfully 
describes net-proton rapidity distributions and their energy and mass dependence. The 
remarkable feature of geometric scaling predicted by the CGC is reflected in the net-baryon 
rapidity distribution, providing a direct test of saturation physics. 

In particular, we show that the peak position in net-proton or net-baryon rapidity dis- 
tributions of centrally colliding heavy ions at ultrarelativistic energies obeys a scaling law 
involving the mass number and the beam energy. We explore the validity of that scaling 
law in comparisons with experimental rapidity distributions from central Pb + Pb and Au 
+ Au collisions at SPS and RHIC energies and find good agreement with the NA49 and 
BRAHMS data on net-proton distributions that were scaled to net baryons. 

We simultaneously investigate net-proton transverse momentum spectra at RHIC en- 
ergies, as measured by the BRAHMS Collaboration, to find reasonable agreement of our 
model results with the data away from midrapidity. There is an increasing discrepancy 
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toward midrapidity, which is expected because the window for saturation effects is reduced 
with decreasing energy and rapidity. 

The parameters saturation-scale exponent A and momentum scale Ql are determined 
from the transverse momentum spectra for Au + Au at 200 GeV and then kept fixed at all 
energies (SPS, RHIC, and LHC) in the calculations for the pt spectrum at 62.4 GeV, and in 
all calculations of rapidity distributions. The third parameter is the overall normalization 
constant C, which we adjust to the data. 

Model calculations for the rapidity distribution in central Pb + Pb collisions at LHC 
energies of 5.52 TeV are obtained; we previously discussed the dependence of the position 
of the fragmentation peak on the gluon saturation- scale exponent A in Ref. [3] and do not 
repeat it here. 

Our analytical scaling law yields an excellent description of the mass dependence of the 
net-baryon distribution at SPS energies in a comparison of S + S and Pb + Pb results. The 
centrality dependence of the Pb + Pb rapidity distribution at 17.3 GeV is well reproduced, 
whereas discrepancies remain for Au + Au at RHIC energies. 

The theoretical result for the mean rapidity loss in ^snn = 200 GeV Au + Au is larger 
than the BRAHMS result as derived from their data, but consistent with the experimental 
upper limit. This emphasizes the importance of a detailed analysis at LHC energies, where 
it may then be possible to determine the value of the saturation-scale exponent A more 
accurately 16) and establish the attainment of gluon saturation. 

To achieve this, forward measurements of identified hadrons — and in particular, 
baryons — for central heavy-ion collisions will be desirable. The ATLAS and Compact Muon 
Solenoid (CMS) detectors for p + p collisions are being extended to the forward region al- 
ready; TOTEM is a dedicated forward detector. It will be useful to exploit these capabilities 
for central Pb + Pb physics as well. 
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